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In this study, we revisit the previous mass relations of mirror nuclei by considering 1/N- and 
1/Z-dependent terms and the shell effect across a shell. The root-mean-squared deviation is 66 keV 
for 116 nuclei with neutron number N > 10, as compared with experimental data compiled in the 
AME2020 database. The predicted mass excesses of 173 proton-rich nuclei, including 98 unknown 
nuclei, are tabulated in the Supplemental Material herein with competitive accuracy. 
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I. INTRODUCTION 


The nuclear mass M(N, Z), where N and Z are the 
neutron and proton numbers of a nucleus, respectively, 
is a fundamental quantity for an atomic nucleus and 
is crucial in cosmology and astrophysics. Several the- 
oretical methods exist for describing and predicting nu- 
clear masses [1-3]. Examples include global models such 
as the Skyrme-Hartree-Fock—Bogoliubov theory [4-6], 
the finite-range droplet model [7-9], the Weizsäcker- 
Skyrme mass formula [10-13], and the relativistic con- 
tinuum Hartree—Bogoliubov theory [14], as well as lo- 
cal models such as the Audi-Wasptra extrapolation 
method [15-19], the Garvey—Kelson mass relations [20— 
24] and their improvements [25-29], and mass relations 
based on neutron—proton interactions [30-33]. Further- 
more, machine learning is widely used to study nuclear 
masses [34-38] and other physical quantities [39-42]. 

Additionally, researchers have developed some formu- 
las based on isospin symmetry to predict the masses of 
mirror nuclei, such as the isobaric multiplet mass equa- 
tion [43-46] and the improved Kelson-Garvey mass re- 
lations [47, 48]. Recently, a simple relationship corre- 
lating the Coulomb energy and shell effect was derived 
to describe the masses of mirror nuclei with a root- 
mean-squared deviation (RMSD) of approximately 200 
keV [49]. This relationship was further improved in other 
studies [50-54] and has expanded to unbound systems 
beyond the proton drip line [53]. 

The purpose of this study is to improve the mass re- 
lations of mirror nuclei presented in Ref. [53] by inves- 
tigating the 1/N- and 1/Z-dependent terms and a sim- 
ple shell correction. The improved relations are used to 
describe the experimental mass excesses of proton-rich 
nuclei with N > 10 with remarkable accuracy and to 
predict some unknowns. The remainder of this paper is 
organized as follows: In Sect. II, we explain the improve- 
ments proposed using empirical formulas for one-nucleon 
separation energies. In Sect. III, we investigate the pre- 
dictive power of the improved mass relations for mirror 
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nuclei. Finally, we conclude this paper in Sect. IV. 


II. IMPROVED MASS RELATIONS OF MIRROR 
NUCLEI 


We begin our discussion with the mass relations of 
mirror nuclei presented in Ref. [53], which are defined as 
follows: 

Ai(N —k, Z) = [M(N —k-1,Z) — M(N —k, Z)] 
—|[M(N,Z—k-—1)—- M(N,Z—k)] 
= S,(N — k, Z) — S,(N,Z—k) — AEnp 
=a + (—)Y F(A- k) t -C , 


c 


(1) 


Ap(N — k, Z) = [M(N — k,Z — 1)— M(N — k, Z)] 
[M(N —1,Z — k) — M(N, Z — k)] 
= Sp(N — k, Z) — Sa(N, Z — k) + AEnp 
= að — (—) B(A- k)! +C 


(2) 
for N = Z, where k is an integer in the range of 1-4, 
Sn (Sp) is the one-neutron (one-proton) separation en- 
ergy, AEnp ~ 0.782 MeV represents the atomic mass 
difference between a neutron and proton, A — k is the 
mass number (A = N + Z), and C is an adjustable con- 
stant. The terms with the parameter a. correspond to 
sophisticated Coulomb energy terms [51], where 


5! =(k+1)(A—k 


c 


2)(A — k = 1)71/ 


(3) 


k(A — k— 1)(A— k)’ — 0.808 , 
= (k—1)(4— k — 2)(A — k — 1 (4) 
k(A —k—1)(A— k)~*3 + 0.808 , 


and 6 in Eqs. (1)-(2) is an optimized parameter that 
represents empirical odd-even staggering in the Coulomb 
energy [52]. 

Deviations between the theoretical and experimental 
values of An [Ap] calculated using Eq. (1) [Eq. (2)] based 
on the AME2020 database [19] are shown in Fig. 1 (a) 
[(b)]. The corresponding RMSDs and number of pairs of 
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Fig. 1. (Color online) Panels (a) and (b) correspond to devia- 
tions between theoretical and experimental values of Ay, and 
Ap calculated using Eqs. (1)—(2) [53], based on AME2020 
database [19]. Red circles represent cases with Z > 20 and 
N — k < 21 (or 20) for An(N — k, Z) [or Ap(N — k, Z)], and 
blue triangles represent cases with Z > 28 and N — k < 29 
(or 28) for An(N — k, Z) [or Ap(N — k, Z)]. 


TABLE 1. RMSDs (in keV) of A, and A, with N —k > 10 
and 1 < k < 4, based on AME2020 database [19]. Second col- 
umn corresponds to improved mass relations of mirror nuclei 
presented by Eqs. (16)—(17), and third column corresponds 
to Eqs. (1)—(2) [53]. Last column M denotes number of pairs 
of mirror nuclei considered. Last row (labeled by A) shows 
results obtained using unified parameters of An and Ap. 


This study Eqs. (1)-(2) N 
An 59 70 46 
Ap 63 80 70 
A 66 78 116 


mirror nuclei considered (denoted by M) are presented 
in the third and last columns of Table 1. Here, five 
experimental data with uncertainties greater than 100 
keV are excluded in our calculations. 

To improve the accuracies of Eqs. (1)-(2), we consider 
empirical formulas for the one-neutron separation energy 
Sn and one-proton separation energy Sp [55, 56], i.e., 


Z 
Sh = (a + aA?) N + a3 + Opair + Osh ; (5) 


N 
Sp = (a + as") z + ag + Opair + Ôsh + Ôcoul ; (6) 


where coul = aeZA~'/3 is the Coulomb term; dpair = 
+a,A7!/? [the signs “+” and “—” correspond to Sn (Sp) 


of a nucleus with an even N (Z) and odd N (Z), respec- 
tively] is the pairing term; dsh = asn is the empirical 
correction of the shell effect with n equaling 0, 1, 2, 3, 
and 4 for Sa (Sp) of a nucleus with N (Z) in the ranges 
of 1-28, 29-50, 51-82, 83-126, and above 127, respec- 
tively; and a; are optimized parameters. The authors 
of [55] highlighted that Eqs. (5)-(6) show an N/Z de- 
pendence similar to the equations for Sp and Sp derived 
from the approximate major-shell lowest-seniority mass 
equation [57], which is based on shell-structure calcula- 
tions. 

Based on Eqs. (1)—(2), A, and A, are related to the 
difference between Sn and Sp. Therefore, in reference 
to Eqs. (5)-(6), the terms dependent on Z/N and N/Z 
should remain, in addition to the Coulomb-energy and 
constant terms, which are well considered in A, and Ap. 
Hence, the correction terms are written as 


Z 
kn _ 
a — N-k’ 
k Z 1/3 ) 
no" (Ack 
yn = —“_(A-h) 
for An, and 
N-k 
iP = , 
N-k (8) 


ô = = (A-k) 


for Ap, with N = Z. In Fig. 2 (a) [(c)], the experimental 
mass differences A,, (Ap) are plotted vs. 6{" (0P) based 
on the AME2020 database [19], where the black squares, 
red circles, blue triangles, and green stars correspond to 
the results of k = 1, 2, 3, and 4, respectively. Clearly, 
the growth rates of A, and A, differ depending on the 
value of k. 

To avoid this difference in k, which results in excessive 
parameters, we consider the neutron—proton interaction 
between the last neutron and last proton (denoted by 
dVin-1p) [32] of the two-mirror nuclei recommended in 
Ref. [50]. Here, ôVin—1p is expressed as 


bVin-1p(N, Z) = -—M(N, Z) — M(N -1,Z - 1) 
+M(N-1,Z)+M(N,Z-1). 
(9) 
Suppose that ôVın—-ıp of the two-mirror nuclei are 
equal [50, 58]; therefore, we have 6Vip—ip(N — k, Z) > 
dVin-1p(N, Z — k). This equation can be rewritten as 


[M(N — k — 1, Z) — M(N — k, Z)] 

— [M(N,Z —k-1)- M(N,Z-k)] ~ (10) 
[M(N—-k-1,Z-1)-M(N-k,Z-1)| 
M(N -1,Z-—k-1)-M(N-1,Z-k)] 
M(N —k,Z — 1)— M(N — k, Z)] 
[M(N -1,Z —k)- M(N,Z —- k) ~ a 
M(N —k-1,Z-1)-M(N-k-1,Z)] 
[M(N-1,Z-k-1)—-M(N,Z—-k-1)] . 
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Fig. 2. (Color online) Panels (a) and (b) correspond to experi- 
mental mass differences An vs. ôf” and 6?; and panels (c) and 
(d) correspond to experimental A, vs. ôP and ôP, respec- 
tively, based on AME2020 database [19]. Black squares, red 
circles, blue triangles, and green stars correspond to k = 1, 
2, 3, and 4, respectively. Different growth rates of A, and 
Ap shown in panels (a) and (c) are unified to a compact tra- 
jectory, as plotted in panels (b) and (d), respectively. 


Based on the definitions of An and A, presented in 
Eqs. (1)—(2), Eqs. (10)—(11) are equivalent to 

Ai(N —k,Z) > 

A, (N —k, Z) ~ 


A,(N —k,Z—1), 
(12) 
A,(N—-—k-1,Z), 
which implies that A, is approximately independent of 
Z for N —k, and A, is approximately independent. of 
N — k for Z, when k = 1-4. Thus, the correction terms 
in Eqs. (7)-(8) can be rewritten as 


T 1 
oa N: k’ (13) 
for Ay, and 
1 
ee 
Oy 7, 7 , (14) 
= zA- k) 


for Ap. Based on the AME2020 database [19], the ex- 
perimental mass differences A, (Ap) are plotted vs. 6? 
(ô?) in Fig. 2 (b) [(d)]. The growth rates of A, and Ap 
are unified into a compact trajectory, as shown in pan- 
els (b) and (d), and are independent of k. This result 
is similar to the finding that strong linear correlations 
between A, (^p) and 62 (6?) are independent of k, as 
mentioned in Ref. [50]. Thus, A, (Ap) with different k 
can be calculated using unified parameters. 


The other correction in this study pertains to the shell 
effect across a shell. Similar to the empirical correction 
of the shell effect for separation energies presented in 
Eqs. (5)-(6), we perform a simple correction for nuclei 
with neutron and proton numbers in different shells, i.e., 


ashi, Z > 20 and N — k < 21 (or 20) for 
A,(N — k, Z) [or A (N — k, Z)| 
ch = 4 Gsh2, Z > 28 and N —k < 29 (or 28) for , 
A, (N — k, Z) [or A(N — k, Z)] 
0, others 


(15) 


where asnı and agp2 are adjustable constants. Owing to 
the limitations of the experimental data, our calculation 
is restricted to nuclei with Z < 50. Thus, only the pa- 
rameters for 20 and 28 shells are presented in Eq. (15). 
Considering the two types of corrections discussed 
above [Eqs. (13)—(14) and Eq. (15)], our improved mass 
relations for mirror nuclei are written as follows: 


Ay(N — k, Z) = ac? + (-)" "(A — ke) 


16 

+ að? + a205 + 6g, — C ( ) 

Ay(N ~~ k, Z) = acoP — (—)76(A = ky (17) 
að? a205 a } C, 


for N = Z, where a;i, a;, 8, and C are optimized param- 
eters; ô and ô are presented in Eq. (3) and Eq. (4), 
respectively. 

Based on the AME2020 database [19], the RMSDs of 
our improved An(N —k, Z) and A (N —k, Z) [Eqs. (16)- 
(17)] with N—k > 10 are presented in the second column 
of Table 1. To reduce the total number of parameters, 
we assume the same parameters for A, and Ap, and the 
results are presented in the last row of Table 1, labeled by 
A. Unification reduced the total number of parameters 
from 14 to 7, with a slight change in the RMSD. The 
results obtained using Eqs. (1)-(2) [53] are provided for 
comparison. As shown, our improvements were highly 
efficient, with the RMSD decreasing from 78 to 66 keV 
for 116 pairs of mirror nuclei. 


Table 2 presents the optimized parameters of A ob- 
tained via least-squares fitting. As shown, the value 
of parameter a, is reasonably close to its typical value 
of 0.70-0.72 MeV, whereas the value of parameter C is 
much smaller than the expected value of 0.782 MeV be- 
cause its contribution is included partly in our correction 
terms. Notably, the signs of parameters asnı and asnz are 
opposite. This is consistent with the results presented in 
Fig. 1, where the red circles (representing cases of asnı) 
are almost above (below) the dashed gray line for A, 
(^p), and the blue triangles (representing cases of asn2) 
indicate the opposite. 


TABLE 2. Unified parameters (in keV) of improved A,(N — 
k, Z) and A,(N — k, Z) [Eqs. (16)—(17)] for nuclei with N — 
k > 10 and 1 < k < 4, based on AME2020 database [19]. 
ac B aı a2 ashi ash2 C 
678 2198 537 —1761 —122 43 317 


Ill. PREDICTION OF MASSES FOR PROTON-RICH 
NUCLEI 


Our improved mass relations for mirror nuclei allow 
one to predict unknown masses of proton-rich nuclei with 
10 < Z < 50. Based on Eqs. (1)—(2), we obtain 
M*(N —k—1,Z) = M(N—k,Z)+ M(N,Z—k-—1) 

— M(N,Z — k) + At? (N — k, Z) , 


(18) 

M"(N —k,Z) = M(N —k,Z—1)+M(N,Z—k) 
M(N -1,Z—k) — A(N — k, Z), 
(19) 


with 1 < k < 4, where AtP(N — k, Z) and Ath(N — k, Z) 
are calculated using Eqs. (16)-(17) with the unified 
parameters. 

Based on Eqs. (18)-(19), up to two masses (M!") can 
be predicted for a nucleus, and their average value with 
the weight of the theoretical uncertainties (o;) is consid- 
ered as the predicted mass M*}, i.e., 


The theoretical uncertainty of Mt}? is defined as o = 
VF [50]. Here, (o;)? (i = 1,2) is the squared RMSD of 
A plus the sum of squared experimental uncertainties of 
the masses of the nuclei involved in the calculation. 

To investigate the predictive powers of Eqs. (18)— 
(20), we predicted the masses of proton-rich nuclei based 
on the AME1995 [15], AME2003 [16], AME2012 [17], 
and AME2016 [18] databases, and then compared the 
predictions with experimental values obtained from the 
AME2020 database [19]. The RMSDs (denoted by 095, 
003, 012, and c16, respectively) and the corresponding 
number of predictions considered (denoted by Npre) are 
presented in the third and last columns of Table 3, 
respectively. The results in (out) parentheses corre- 
spond to cases that include (exclude) experimental data 
with uncertainties greater than 100 keV. The RMSDs of 
Ref. [53] are listed in the second column for comparison, 
which indicate the competitiveness of our predictions. 
Based on the AME2020 database [19], the mass excesses 
of 173 proton-rich nuclei (including 98 unknowns) with 
10 < Z < 50 were predicted and are tabulated in the 


TABLE 3. RMSDs (in keV) of predicted masses of proton- 
rich nuclei compared with experimental data obtained from 
AME2020 database [19]. 095, 03, 12, and oig correspond 
to extrapolation based on AME1995 [15], AME2003 [16], 
AME2012 [17], and AME2016 [18] databases, respectively. 
Last column Nore denotes number of predictions consid- 
ered. Results in (out) parentheses correspond to cases includ- 
ing (excluding) experimental data with uncertainties greater 
than 100 keV. Results of Ref. [53] are listed in second column 
for comparison. 


Ref. [53] This study Nore 
T95 113 (146) 90 (133) 19 (25) 
503 116 (210) 95 (198) 20 (20) 
O12 76 (72 77 (68) 11 (9) 
716 84 (79 84 (74) 8 (7) 
E O AME2020 O Thiswork A Ref. [53] @ Ref. [62] 
0.4| (a) 8s (b) Co (c) *Cu (d) "Kr 
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Fig. 3. (Color online) Nuclear masses deviating from those 
in AME2020 database[19]. Solid (hollow) black squares cor- 
respond to experimental (estimated) values obtained from 
AME2020 database [19]. Masses predicted in this study and 
Ref. [53] are labelled by hollow red circles and blue trian- 
gles, respectively, and recent measurements obtained from 
Ref. [62] are labeled by solid green diamonds. Experimental 
and theoretical uncertainties are labeled as well. 


Supplementary Material [59] along with the correspond- 
ing theoretical uncertainties. 

In Fig. 3, the predicted masses of eight nuclei are pre- 
sented (labeled with red circles) with respect to those 
presented in the AME2020 database [19]. Panels (a)—(e) 
correspond to five nuclei with experimental uncertainties 
greater than 100 keV [19], whose experimental masses 
are excluded from our calculations, as mentioned in Sec- 
tion II; panel (f) corresponds to the nucleus, in which the 
unbound nature is confirmed experimentally [60, 61]; 
and panels (g) and (h) correspond to nuclei with recent 
new measurements (labeled by green diamonds) [62]. 
The predicted masses obtained from Ref. [53] (labeled 
with blue triangles) are provided for comparison. Our 
predicted masses are more similar to the experimental 
values and show significantly lower theoretical uncertain- 
ties, except for °°Cu in panel (c). 

For ĉ!Ga in Fig. 3 (g), the recent experimental mass is 
—47114 (12) keV [62], which is three times more precise 


TABLE 4. Newly measured nuclei masses (in keV) obtained from Ref. [63] and their predictions obtained from Ref. [52], 
Ref. [53], Ref. [54], and our improved relations [Eqs. (18)—(20)]. Values in parentheses indicate corresponding experimental 
or theoretical uncertainties. o{ and gô (in keV) in last two rows correspond to RMSDs between theoretical and measured 
values [63] for all 12 nuclei and the same six nuclei presented in Ref. [52], respectively. 


Nuclei Ref. [63] Ref. [52] Ref. [53] Ref. [54] This work 
5872n —42248(36 —42319(82 —42334(68) —42327(41) —42327(47) 
Ga —40034(46 —39995(84 —40030(71) —39982(41) —40022(50) 
Ga —47168(21 ~ —47099(97) —47085(44) —47129(68) 
Ge —42289(37 —42332(86 —42328(73) —42349(41) —42373(54) 
Ge —46978(15 ~ —47010(96) —46993(41) —47022(66) 
S4As —39710(110) —39662(93 —39512(82) —39562(41) —39601(62) 
65 As —46806(42 - —46751(96) —46760(41) —46790(66) 
6 Se —41982(61 —42046(98 —41847(88) —41890(41) —41985(71) 
87Se —46549(20 - —46581(96) —46588(41) —46604(66) 
Kr —41320(140) —41560(101) —41269(96) —41333(41) —41320(78) 
Ky —46056(24 - —46037(96) —46068(42) —46070(66) 
Sr —46200(150) - —46302(96) —46356(41) —46346(66) 

ot = - 86 81 67 

oh = 110 107 85 65 


than the experimental value provided in the AME2020 
database [19], i.e., —47134.7 (38.0) keV. The mass of 
60Ga, as shown in Fig. 3, (h) was recently measured to 
be —40005 (30) keV [62], which is more than 400 keV less 
than the estimated value presented in AME2020 [19], i.e., 
—39590# (200#*) keV. Our predictions for the masses of 
both nuclei are more similar to their new measurements, 
in comparison with the values presented in Ref. [53]. 

Table 4 lists the predicted masses and those presented 
in Ref. [52], Ref. [53], and Ref. [54] for comparison with 
the newly measured experimental values obtained from 
Ref. [63]. The corresponding RMSDs are listed in the 
last two rows of Table 4, where o{ and o4 denote the 
RMSDs for all the 12 nuclei and the same six nuclei in 
Ref. [52], respectively. Our predictions agree well with 
the experimental values, and the RMSD can be further 
reduced to 48 keV for nine nuclei, excluding ® As, “Kr, 
and “°Sr, whose experimental uncertainties are greater 
than 100 keV. 


IV. SUMMARY 


In this study, we revisited the mass relations of mirror 
nuclei proposed in Ref. [49] and improved in Ref. [53] 
by considering 1/N- and 1/Z-dependent terms and the 
shell effect across a shell, which originated from the em- 
pirical formulas of one-nucleon separation energies. The 
improvements of our formulas were shown to be effective 
for nuclei with neutron number N > 10, as indicated 
by an RMSD of only 66 keV for 116 pairs of mirror nu- 
clei, as compared with the experimental data presented 
in the AME2020 database [19] and 90 keV for 19 proton- 
rich nuclei extrapolated from databases AME1995 [15] 
to AME2020 [19] (excluding experimental data with un- 
certainties greater than 100 keV). 


The predicted mass excesses of proton-rich nuclei with 
10 < Z < 50, based on the AME2020 database [19] and 
our improved formulas, as well as the corresponding the- 
oretical uncertainties, were tabulated, as presented in the 
Supplementary Material [59]. These predictions agreed 
well with the newly measured nuclei masses [63], with 
an RMSD of only 67 keV for 12 nuclei. The RMSD was 
further reduced to 48 keV by excluding three experimen- 
tal data, with uncertainties greater than 100 keV. We 
believe that these predictions of proton-rich nuclei pro- 
vide valuable information for investigations pertaining 
to nuclear physics, such as two-proton radioactivity [64]. 
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